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Abstract 

We consider the scalar-tensor theories of gravity extended by the pseudoscalar couplings 
to matter and gauge fields and derive constraints on the CP-odd combinations of scalar and 
pseudoscalar couplings from laboratory spin precession experiments and from the evolution 
of photon polarization over cosmological distances. We show the complimentary character of 
local and cosmological constraints, and derive novel bounds on the pseudoscalar couplings 
to photons from the laboratory experiments. It is also shown that the more accurate 
treatment of the spin content of nuclei used in the spin precession experiments allows to 
tighten bounds on Lorentz-violating backgrounds coupled to the proton spin. 



1 Introduction 



The discovery of dark energy [1] instigated many developments in cosmology and particle 
physics during the last decade. To date all observational data are consistent with the most 
economic possibility: the dark energy is just a cosmological constant, and as such does 
not evolve over the cosmological time scales. On the other hand, it is intriguing to think 
about the alternative explanations related to a drastic change of the infrared physics. In 
parallel to the attempts of modifying gravity on large scales [2] , there is a renewed interest 
in the cosmological scalar fields that are nearly massless, and manifest themselves as a 
"dark energy" component over large cosmological distances [3]. 

An interesting twist to the well-known story of cosmological scalars comes from the pos- 
sibility of their interaction with matter and gauge fields (For purely cosmological signatures 
of "interacting" quintessence, see e.g. [4]). In fact such theories exhibit a rich plethora of 
phenomena that go beyond pure cosmological effects, which we would like to illustrate on 
the following toy example. Let us consider a Lagrangian for the scalar field interacting 
with a Standard Model fermion ip {e.g. electron) and a gauge field ^4^ {e.g. photon), 

C = ^d,<l)dy - \/(0) + ^{tD^-f^ - m)ilj - ^F^uF^-" (1.1) 
-cs^4>ipip - cp^4>ipij5ip - cs'y4>Ffj_^F'"' - cp^(t)F^^F^'' . 

Here csi and Cpi parametrize the strengths of the scalar and pseudoscalar couplings to 
photons and fermions, while F^^, and F^n denote the usual and dual field strengths, and 
Dfj, is the covariant derivative. Written in flat space, Lagrangian (1.1) can be trivially 
generalized to curved backgrounds, and to nonlinear couplings to matter. Starting from 
(1.1), one can immediately infer a number of interesting consequences, a partial list of which 
is given below. 

1. The existence of a new long-range force distinguishable from spin-two gravity. The 
scalar field contributes to the gravitational force, adding ~ on top of the familiar 
Newton constant mediated by gravitons. Such a force leaves distinguishable imprints 
via relativistic corrections and/or composition dependence (effective violation of the 
equivalence principle) . 

2. The existence of a preferred Lorentz frame associated with dt(f). If is a very light 
quintessence-like field, then there is a preferred frame where cosmologically = 
{(f), 0, 0, 0). For most of the models this frame coincides with the frame of the cosmic 
microwave background (CMB), and |0| is limited by {pd.e.{^ -|-w))^/^, where w is the 
dark energy equation of state parameter. 

3. Variation of masses and couplings in time and space. Effective values of masses and 
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coupling constants vary in space and time, mphys{t, x) = m-\-cs'^(t>{ti x), following the 
^-profile. 

4. Coupling of polarization to velocity relative to the CMB frame. A particle moving 
relative to the CMB frame acquires a helicity-dependent interaction, Hint ~ (Sn)0, 
where n is the direction of propagation. This way, the Cp-y-proportional interaction 
would result in the rotation of polarization for photons propagating over varying 
0-background. 

5. Photon-scalar conversion. In the presence of an external electromagnetic field a pho- 
ton can "oscillate" to a quantum of the scalar field thereby e.g. reducing the lumi- 
nosity of distant objects or providing additional channels for star cooling. 

6. Coupling of spin to the local gravitational force. Scalar coupling gsip will lead to the 
local field gradient V0 generated by massive bodies, which is closely parallel to the 
vector of local free-fall acceleration. The pseudoscalar couphngs then create a Zeeman- 
like splitting for the spin of ■0-particles in the direction of the local gravitational 
acceleration. Hint ~ (Sg). 

It is remarkable that such a simple Lagrangian leads to a number of quite different 
phenomena. Unfortunately, at this stage the exciting phenomenology of "interacting dark 
energy" lives in a pure theoretical realm: there is no confirmed experimental evidence for 
any of the effects on our list^. Consequently, there are only upper limits on the combinations 
of the couphngs in Lagrangian (1.1) that can be quoted. Nevertheless, many of the effects 
on our list have found an extensive coverage in the theoretical works. Most notably, the 
changing couplings were discussed, for example, in Refs. [9, 10, 12], the photon-scalar 
conversion considered in Refs. [13], and the fixed frame effects versus the cosmological 
evolution of photon polarization were addressed in a series of papers [14, 15, 16, 7, 8]. For 
the limits on scalar-induced corrections to gravitational interactions we refer the reader to 
recent reviews [17] and references therein. In contrast, the last item on our list, the spin 
coupling to the local gradient of the scalar field received far less attention (see e.g. [18]). 

The purpose of this paper is to show that the pseudoscalar couplings of the Brans- 
Dicke type scalar can indeed be subjected to stringent laboratory constraints that are 
complementary to the cosmological limits. The high-precision spin precession experiments 
constrain pseudoscalar interactions both in the fermion and photon sectors. In the rest of 
this paper we present the set-up for our model, briefiy review the effects created by the 
cosmological evolution of 0(t), investigate the local spin effects created by the gradient of 
0, and set the limits on the admissible size of the pseudoscalar couplings. 

-'^A tantalizing hint on the redshift evolution of the fine structure constant was reported in Ref. [5], which 
so far has not been corroborated by other searches [6] . Also, an earlier claim of the non-zero pseudoscalar- 
induced anisotropy in polarization signal [7] was disputed in the literature [8]. 
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Before we delve into studying the physical effects induced by the pseudoscalar couplings, 
we would like to add a word of caution addressed to all models of " interacting quintessence" . 
The models of light scalar fields represent a formidable challenge at the quantum level, 
as there are no fundamental reasons for a scalar to remain massless or nearly massless. 
The scalar interaction of such field makes the whole problem even more difficult, if not 
impossible, from the point of view of "technical naturalness": the loops of Standard Model 
(SM) fields tend to generate big corrections to V{(j)) even with a relatively small ultraviolet 
cutoff parameter, which would be in conflict with requirements ~ H [10, 15, 19]. There 
is no clear resolution to this problem, which essentially prevents the fully consistent study 
of dynamics. Instead, one has to rely, perhaps too optimistically, that the problem of 
near-masslessness of the scalar held could be cured by the same mechanisms that make 
the cosmological constant small and meanwhile keep V{(j)) fixed by hand. To finish this 
"disclaimer" on an optimistic note, we would like to remark that the pseudoscalar couplings 
do not make this problem worse. Indeed, in essence the pseudoscalar couplings give only 
derivative interactions, and therefore do not affect the potential V{4>) at perturbative level. 



2 Adding spin couplings to scalar-tensor theories 

We would like to formulate our reference Lagrangian at the normalization scale just below 
the QCD scale, so that the effective matter degrees of freedom are electrons, photons, 
nucleons and neutrinos. Splitting the 0-field Lagrangian into the scalar and pseudoscalar 
parts, 

/: = jrs + J^P: (2.2) 
we choose the following parametrization, 

Cs^ld,(l>d^(l>-V{<j>)- E l^rn,^P0, ~ -^F,,F^\ (2.3) 

2 j=e,p,n ^Sj Msj 

and 



>Cp= E ^ ^.7m75^. - ^i^M.^'^^ (2-4) 

j=e,p,n,v 



Mpj ""^^ Mp^ 



Lagrangian (2.4) includes all possible pseudoscalar interaction at mass dimension five level. 
Notice that the pseudoscalar interactions can be chosen in a slightly different form, 'i/^'y^ijj, 
as in (1.1). This does not mean, however, that the set of our operators should be enlarged. 
The two type of operators, pseudoscalar and axial vector, are related on the equations of 
motion. These equations are in general anomalous, but since we include the interaction with 
FF explicitly, we can assert that Lagrangian (2.4) is indeed complete in a given dimension 
of the operators. 
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The scalar part of the Lagrangian (2.3) leads to new contribution to gravitational force, 
and to change of masses and couplings. Since in this paper our main interest is in spin 
effects, we are going to make simplifying assumptions of approximate universality of the 
0-mediated attractive force, 

Mse^Msp^Msn = Ms, and Ms^:^Ms. (2.5) 

At distances shorter than the Compton wavelength of 0-quanta the Newton constant re- 
ceives contributions from both spin-two and spin-zero exchanges, 

G„ = G»„(l + H|l), (2.6) 

where G^^ is the unperturbed Newton constant due to graviton exchange, and Planck mass 
is defined as Mpj = (87rG'^)-V2 = 2.4 x lO^^GeV. 

If needed, the pseudoscalar couplings could be "hfted" from the nucleon level to the level 
of individual quarks. Using the experimental results for the spin content of the nucleon 
combined with SU (3)-flavour relations [20] one gets 

Mpl ~ Q.mpl - 0.4Mp] - O.lMpi; Afp^ ~ 0.8Mp] - 0.4Mp^ - 0.1Mp,\ (2.7) 

where the light quark couplings are normalized at the scale of 1 GeV. 

Using the appropriate field content, one can determine the renormalization group evolu- 
tion of the pseudoscalar couplings. In general the equations governing this evolution takes 
the following form, 

where the logarithm is taken between the ultraviolet scale A and the infrared scale 
Latin indices indicate fermionic fields and Greek indices indicate the gauge bosons of the 
SM group. The renormalization group coefficients aj^, feja, c^j, and d^p depend on charge 
assignments and coupling constants of field running inside the loops. The precise form of 
these coefficients is not of immediate interesting to us, but we would like to emphasize the 
following important observation: at any loop level the derivative couplings to fermions do 
not generate couplings to F^^F^" . In other words, 

Cai = 0. (2.9) 

Whatever size of the pseudoscalar couplings between photons and is generated by some 
(perhaps anomalous) ultraviolet scale physics at energies order A, it is preserved by the 
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subsequent evolution to the lower scales. In fact, this refers both to the logarithmic running 
and to the threshold corrections. This observation dehneates two important classes of 
models: there are models where both fermion and photon pseudoscalar couplings present 
in the Lagrangian, and there are models where only couplings to fermions are present. The 
models where couples only to gauge bosons would necessarily be fine-tuned, as quantum 
effects in (2.8) would definitely generate induced couplings to fermions. 

Existing constraints on the model can be devided into pseudoscalar and scalar con- 
straints. The constraints on the universal scalar coupling Ms can be derived from the 
constraint imposed by the Cassini satellite data on the post-Newtonian parameter 7 [21], 

I7I < 4 X 10"^ =^ Ms> 400 Mpi. (2.10) 

The constraints on the non-universal part of the scalar coupling are several orders of mag- 
nitude stronger. The scalar coupling to photons is constrained via the limits on the time 

variation of the coupling constant, and less directly via the composition-dependent contri- 
bution to local accleration. Typically, one has Ms-y > lO'^Afpi. In contrast, the pseudoscalar 
couplings are far less constrained. The leading source of constraints arc the energy losses 
mechanisms in stars [22] , and for electrons, photons and nucleons all constraints are in the 
ballpark of 

\Mp\ > (10^° - 10^^) GeV ~ (10-^ - 10-^) x Mpi. (2.11) 

In the next section, we are going to show that if both pseudoscalar and scalar couplings are 
present, some constraints on Mp can be significantly improved. 



3 Cosmological constraints on the model 

To derive cosmological constraints on pseudoscalar couplings wc remind the reader that the 
presence of a time-evoloving scalar field with a pseudoscalar coupling to photons leads to a 
rotation of polarization for photons. The resulting angular change in the linear polarization 
for a photon propagating from point 1 to point 2 is simply related to the change of between 
the two points. 

Following the work of Carroll [15] and the original analysis of Ref. [24], we use the limit 
on the extra rotation of polarization from distant source (3C 9) at redshift z = 2.012 as 
|A^| < 6°, 

|0(z = 2)-0(z = O)| 



Mp^ 



< 0.052. (3.2) 



Even more distant sources of polarization are available in the studies of the cosmic 
microwave background. The £'-mode polarization map of the sky has been produced [23] , 
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which agrees well with the expectation based on the temperature map. This constrains 
the amount of extra rotation of polarization introduced by (f)FF interaction between the 
surface of last scattering and z = 0. Recent numerical analyses of the CMB data provide a 
constraint on the amount of extra rotation at the level of |A^^| < 6° [16] (the same limit of 6° 
is purely coincidental), which allows to extend (3.2) to the redshifts of photon decoupling, 

Zdec — 1100, 

|0(..ec)- 0(^ = 0)1 ^ (3.3) 

Finally, we would hke to point out that the CMB polarization signal is generated in the 
narrow window of redshifts that correspond to the "last scattering" surface, and therefore 
the existing measurements constrain the amount of extra rotation within the thickness of 
this surface, 

^ |A0(z,ec±AW2)| <0(1), (3.4) 



where Az^ec — 200 correspond to the thickness of the last scattering surface. The violation 
of this bound would suppress the strength of polarization signal, which is well measured. 

With these bounds at hand, we are ready to translate them into the constraints on the 
parameters of our model. However, the cosmological constraints depend very sensitively 
on what we assume about the scalar couplings of to dark matter and even more so on 
the choice of the potential T^(0). Since the number of options is infinite, we would like to 
consider in detail two well-motivated cases. 

Case 1. The simplest case is when the potential for is nearly flat and the evolution of (p 
is slow. In this case one can linearize V{4>), 

+ (3.5) 

where pA is approximately equal to the measured value of dark energy density, and M\ is 
a new parameter on the order of the Planck scale and/or Ms- In the limit when the back- 
reaction of on Friedmann's equations is neglected one can find an analytic expression 
for the evolution of in the flat Universe [10]. In this approximation the time evolution 
of the scale factor can be expressed via the scale factor and the Hubble parameter today 
{tnow = ^o): Hq — H{t — to) — d/a\t=to and a{t — to) = oq, as well as the current energy 
densities of matter and cosmological constant relative to the critical density, = Pm/ Pc 
and flA ^ Pa/ Pc- 



(3.6) 



The equation of motion for the scalar field receives forcing terms directly related to dark 
energy and matter densities: 



^ ^ Ms Ma 



Ms \ aj Ma 



(3.7) 
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where we made an assumption of the universal strength of (f) couphng to matter, including 
the dark matter. This equation can be integrated out explicitly [10] to give 



\i ' 


--) 


_V2Ma 


Ms) 



sinh(6i) 



' ^ " Ms sinh(6io) 

where the following notation has been introduced: 



(3.8) 



b = ^nlHo. (3.9) 

This solution implies boundary conditions (j)\t^o ~ not too large and 0(to) = 0. The 
first condition is automatically satisfied as does not evolve rapidly during the radiation 
domination, and the second condition is simply a choice, possible since (f) enters linearly in 
the Lagrangian. It is easy to see that in the hmit oit <^to the dependence of cf) on redshift 
is logarithmic, 

(f){a) ~ const — ^ ln(a/ao) at teq <^t <^ to. (3.10) 

Ms 

Now we can use the evolution (3.8) to impose limits on the combination of Mp and Ms{a) 
parameters using observational constraint (3.2). We do it for three separate representative 
cases: for the equal couplings to the matter and dark energy density, and for couplings to 
dark energy and matter only: 

36 M|i; (3.11) 
30 M^j; (3.12) 
6.1 M^i. (3.13) 

These limits generalize the analysis of Ref. [15] where only the F(0)-induced optical rotation 
was considered. We also notice that both (3.11) and (3.12) are about one order of magnitude 
stronger than (3.13), which is the consequence of (ao/a)^ ~ 8 enhancement of matter density 
over the cosmological constant at redshifts ~ 2. 

Due to logarithmic dependence on redshifts at t ^ (3.10) there is about one order of 
magnitude gain in the strength of the constraint when using the CMB limit (3.3) for the 
case of finite Ms, 

Ma^ Ms or |Ma| ^ oo =^ \Mp^Ms\ > 255 M^j. (3.14) 

In order to see the maximal sensitivity to Mp, we can saturate the constraint on Ms 
(2.10), which results in Mpj ^ O(Mpi) at maximally allowed Ms- 

Case 2. Going away from the linearized case, we consider the cosmological evolution of 
0-field approaching some local minimum of V{(f)), 

2 

V{<P)^PA + '^{<l>-<l>of. (3.15) 



Ma = 


Ms 


\MpyMs\ 


> 


|Ma| - 


^ oo 


\Mp^Ms\ 


> 


\Ms\ - 


->■ oo 


\Mp^Ma\ 


> 
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If the mass of the field is well above the current Hubble parameter, ^ Ho, then the 
evolution of (p starts long before the present epoch. A well-known solution for (p in this case 
are the oscillations around the minimum with the amplitude that red-shifts as a~^/^. If 
the initial deviation of from equilibrium was (pin at the time tin when oscillations began, 
H{tin) ~ m^, then the subsequent evolution in the radiation domination will be given by 

(pit) - (po + (pinly^j cos[m^{t - tin) + a], (3.16) 

where a is some phase factor. Because of the red-shifted amplitude in (3.16), the constraints 
provided by the CMB are clearly more advantageous than the low z constraints. However, 
the oscillations of (p (3.16) make it difficult to define (p{zdec), and consequently the analyses 
of [16] with limits (3.3) are not directly applicable and instead one should resort to limits 
(3.4). Still, if the initial deviation of (p field from its minimum is on the order or less 
than the pscudoscalar coupling Mp^, and oscillations begin earlier than the decoupling, 
then the cosmological evolution of polarization provides no constraints on the size of the 
pscudoscalar coupling, 

< |Mp^|; tin^tdec no constraints on Mp^. (3-17) 

This is an important observation, since the first condition \(pin\ < \Mp^\ is quite natural if 
(p field has a phase-like origin similar to e.g. QCD axion, and tin <S tdec is satisfied for all 
masses of (p in excess of 10~^® eV. 



4 Local spin precession constraints 

As we have shown in the two previous sections, the cosmological constraints on pscudoscalar 
couplings apply only to Mp^, and not to fcrmionic couplings. Moreover, all cosmological 
constraints will be eliminated if the field starts oscillating much earlier than the decoupling 
of the CMB photons (3.17). This leaves a large domain of parameter space where only 
the local experiments are going to be sensitive to the pscudoscalar couplings. We wish to 
consider them in this section. Before we do that, we would like to note that the couplings of 
spins to the local gravitational (spin-two) field has been extensively studied in the literature 
[25, 26, 27]. Of main interest for us is the conclusion reached in these works that g • S 
coupling does not arise in general relativity. Therefore, if detected, it can be thought of as 
a distinct signature of the scalar exchange. 

Since most of the experiments deal with non-relativistic atoms and nuclei, it is conve- 
nient to use the non-relativistic Hamiltonian, 
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where a — '^l\S\ — 2S. The local gradient of is one-to-one related to the gravitational 
acceleration, 

V^ = g^. (4.19) 

SO that the strength of interaction of each spin to the gravitational field is given by (yf x 
2M'^y/ [MsMpj). Gravitational acceleration has dimension of energy in particle physics 
units of c = ^ = 1, and corresponds to the frequency splitting of spin up and spin down 
states i/g = 2 X 9.8 10^cm/s^/(27r x 3 10^°cm/s) =10.4 nHz. Unhke most problems in 
quantum mechanics where "up" and "down" are usually a matter of convention, in this 
theory these words should be used literally. Only a handful of spin precession experiments 
ever reached the sensitivity lower than 10 nHz, among them the experiments searching 
for the permanent electric dipole moments of diamagnetic atoms [28], where the statistical 
sensitivity is comparable or better than 10 nHz. Unfortunately, this sensitivity is related to 
the energy difference of spins in parallel and anti-parallel electric and magnetic fields and 
does not translate into the limits on spin interaction with the vertical direction. 

Dedicated search for g • S interaction was pursued in [29] (and earlier in [30]), where a 
~ //Hz accuracy was achieved. In particular, experiment [29] compared the precession fre- 
quencies of two mercury isotope spins, ^^^Hg and ^°^Hg for different orientations of magnetic 
field and set a limit of 2.2xl0~^° GeV for the spin-dependent component of gravitational 
energy. Another group of measurements that can be used to limit the pseudoscalar cou- 
plings are the spin precession experiments that searched for the effects of Lorentz violation 
[31, 32] and the experiment with spin-polarized pendulum [33]. The absence of sidereal 
modulation of spin precession, confirmed by these experiments, sets the limit on the cou- 
pling of spins to any direction in space that does not change as the Earth rotates around its 
axis. Besides useful limits on Lorentz- violating theories [34], such effects will constrain the 
pseudoscalar couplings in combination with V0 created by astronomical bodies other than 
the Earth. The solar contribution to V0 is smaller than Vcf) Earth by a factor of ~ 6 x 10~^, 
thereby reducing the strength of the constraints extracted from sidereal variations by the 
same amount. Putting different results together, and assuming that the range of the force is 
comparable to or larger than the solar system, we arrive at the following set of constraints, 

\Mpr,Ms\ > 5 X 10-^ M|i Ref. [29] (4.20) 
\Mpr,Ms\ > lxlO-^M|i Ref. [31,32] (4.21) 
\MpeMs\ > 2 X 10~^ M^i Ref. [33] (4.22) 

Bounds (4.20) and (4.21) are derived in the assumption of Ref. [35] that the spin of the 
nucleus is given by the angular momentum of the of the outside nucleon, which happens to 
be a neutron for all nuclei used in the most sensitive searches (^He, ^^^Xe, ^^^Hg, ^°^Hg). 
Consequently, the limits are formulated on the pseudoscalar coupling to neutrons, as it is 
also the case for the limits on the external Lorentz- violating axial- vector backgrounds [35] . 
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nucleus 




J, I 










-2.13 


1/2, 


1 


1.04 


-0.04 




-0.78 


1/2,0 


1 


0.76 


0.24 


199 Hg 


0.50 


1/2, 1 


-1/3 


-0.31 


-0.03 


20lHg 


-0.56 


3/2, 1 


1 


0.71 


0.29 



Table 1: Composition of the nuclear spin 



In fact, one can refine these bounds and impose separate constraints on the strength of 
the pseudoscalar coupling for protons and neutrons. Although most of the nuclei in atoms 
used in experiments [29]-[32] have a valence neutron outside of closed shells, one can use 
the information on the magnetic moments of these nuclei together with simple theoretical 
model of nuclear structure to deduce the proton contribution to the total nuclear spin. To 
be specific we shall assume that the magnetic moment of the nucleus is composed entirely 
from the spin magnetic moment of the valence neutron and spin magnetism of polarized 
nuclear core, 

/x = ^„(aW)+/.,(a(^)) (4.23) 

In these equations, /x, /Xp, /x„ are the magentic moment of the nucleus, proton and neutron. 
Numerical estimates show that the orbital contribution to the magnetic moment // in the 
nuclei of interest is less important than the spin contribution since the neutron orbital 
contribution is zero and the proton orbital contribution is small in comparison with the 
proton spin contribution. The latter is enhanced by the large value of the proton magnetic 
moment /ip = 2.8, which justifies the neglection of proton orbital magnetism for low I 
obitals. Neglection of spin-orbit interaction makes total spin conserved and its total value 
equal to the average spin of the neutron above the unpolarized core, {cr^p). The latter is 
equal to 1 for j = / -|- 1/2 and —j/ {j + 1) for j = I — 1/2, where j is the value of the nuclear 
angular momentum, and I is the orbital quantum number of the valence neutron. Using 
these simple formulae (4.23), we determine (cr^"-*) and {a^f^) for observationally relevant 
cases of ^^^Hg, ^°^Hg, ^^^Xe, and ^He as shown in Table 1. 

One can see that the constribution of the proton spin into the total spin of these nuclei, 
especially ^^^Xe and ^°^Hg, can be as high as 30%, and therefore the proton pseudoscalar 
coupling is also limited in these experiments. For example, experiment [29] limits the 
following combination of the proton and neutron couplings: 

\Mp^sMs\ > 1.5 X 10"^ M|i, where Mpjg = -0.5Mp^ + 0.7Mp^ (4.24) 
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The relative enhancement of the proton contribution is due to a rather close cancellation 
of neutron contribution to the differential frequency of spin precession for ^^^Hg and ^°^Hg. 

As a bi-product of our analysis, we can improve the bounds on the Lorentz-violating 
axial- vector couplings in the CoUaday-Kostelecky parametrization [34]. Indeed, the spatial 

components of the axial vector background to protons, 6^, is constrained in the same exper- 
iments, Refs. [31, 32], in particular because of the substantial contribution of proton spin 
to the spin of ^^^Xe. For example, the interpretation of the null result of the most sensitive 
experiment [32] with the use of the analysis [35] that assumes (cr^")) = (cr^°-*), (c^^-*), and 
our work differ in the following way: 

2t:vlv = 265") {i - ^ ^ -3.56S"\ Ref.[35] 
2t^vlv = 2(0.766f ) + 0.2461^^ - ^(1.046^^ - QMhf^)) = -4.26^^ + Q.lllf'^ this work. 



e 



where vlv = 53±45 nHz is the experimentally measured (and consistent with zero) Lorentz- 
violating frequency shift [32]. Obviously, the contribution of llf'' to vlv is non-negligible, 
and imphes that \})f^ \ < few x 0(10"^^) GeV, which is far better than the results of dedicated 
searches of Lorentz violation in the proton sector with e.g. hydrogen maser [36]. 

Besides the constraints on nucleon and electron couplings, the same clock comparison 
experiments allow to set limits on Mp^. For example, for an atom (or nucleus) with the 
total angular momentum J, the matrix element of interaction is not zero. 



where /t is a dimensionless matrix element that can be calculated explicitly. For the ground 
state of the hydrogen atom, k, is given by 

.= ^=1^, (4.26) 

where Oq and /i^ are Bohr radius and magneton, and a is the fine structure constant. 
This calculation take into account the magnetic field generated by the electron magnetic 
moment, and the electric field of the proton. If we consider both E and B created by the 
electron, we discover that the result has a logarithmic divergence in the ultraviolet regime 
that has the interpretation of 1/Mpe being generated by 1/Mp^. Even with a modestly low 
choice of the cutoff, the coefficient is going to be on the order of ck/tt ~ 0(10^^) and thus 
parametrically larger than (4.26). 

What happens if instead of an atomic electron we consider a nucleus where the electric 
field is considerably stronger? To understand the scahng of the effect with Z, we consider 
a simplified case of a single s-wave neutron above the closed nuclear shells with "uniform" 
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distribution of its wave function inside the nucleus, which also has uniform charge distri- 
bution within a sphere of radius Rn — 1.2 fm (A)^/^. The resulting k, can be expressed in 
terms of the neutron magnetic moment, 

K = 5 ^ ^ 1 ^ ^ 0.05 - 0.07 for Z ^ 80, (4.27) 
5 Rn 5 rUpRN 

where the overall numerical coefficient follows from the approximation of the radial matrix 
element, {ry{2R%) - 3/2) 

r<iijv = —6/5. Although an overall numerical coefficient in 
estimate (4.27) cannot be taken very seriously, the parametric dependence on Z, and 
Rn is certainly expected to hold for large nuclei. For mercury this effect is larger than 
the loop-induced admixture of the photon coupling into the nucleon coupling. Thus we 
can deduce the sensitivity of spin precession experiments to the pseudoscalar couplings to 
photons at 5% level from the coupling to neutrons: 

\Mp^Ms\ > 0(10"^) M|i Ref. [29] (4.28) 

One can sec that the combined bounds from the clock comparison experiments arc compa- 
rable to or better than the product of separate bounds (2.10) and (2.11). Unfortunately, 
these bounds do not allow to probe the pseudoscalar coupling to fermions all the way to 
the "natural" scale Mp Mpi. 



5 Conclusions 



Our paper considers the constraints on the combination of scalar and pseudoscalar cou- 
plings in the scalar-tensor theories of gravity. The strongest constraints come from the 
considerations of cosmological evolution of polarized light, and in the best case scenario of 
the maximal scalar coupling, consistent with constraints on Brans-Dicke theories, the sensi- 
tivity to the pseudoscalar coupling to photons can be as large as the Planck scale. However, 
the cosmological constraints are not sensitive to the derivative pseudoscalar couplings to 
fermions as they do not induce corresponding photon couplings even at the loop level. We 
also point out that for a wide range of pseudoscalar masses, one can avoid cosmological 
constraints due to the red-shifting of 0-oscillations. Therefore, the laboratory constraints 
on spin precession from locally generated gradient of (f) are complimentary to cosmological 
bounds. We revisited lab bounds to find that the most sensitive experiments are still few 
orders of magnitude below the sensitivity to Planck-scale-suppressed couplings. We also 
note that the local spin precession experiments provide sensitivity to the pseudoscalar cou- 
pling to photons, through the relatively large matrix element of 0B • E interaction inside 
atomic nuclei. As a separate remark, we have shown that nuclei of atoms used in the high- 
precision clock comparison experiments have significant proton contribution to their spins. 
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This allows to set separate constraints on pseudoscalar couplnigs to neutrons and protons, 
and improve the limit on Lorentz-violating axial-vector backgrounds in the proton sector. 
Further progress in experiments searching for a preferred Lorentz frame would also provide 
better sensitivity to the scalar-tensor theories extended by pseudoscalar couplings. 
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